Pressure dependence of two-level systems in disordered atomic chain 



O 

(N 

< 

oo 



CO 



c 

o 
o 



> 

m 
o 



X 



A. ShelkaiE and V. Hizhnyakov 
Institute of Physics, University of Tartu, Riia 142, 51014 Tartu, Estonia 

The dependence of two-level systems in disordered atomic chain on pressure, both positive and 
negative was studied numerically. The disorder was produced through the use of interatomic pair 
potentials having more than one energy minimum. It was found that there exists a correlation 
between the energy separation of the minima of two-level systems A and the variation of this 
separation with pressure. The correlation may have either positive or negative sign, implying that 
the asymmetry of two-level systems may in average increase or decrease with pressure depending on 
the interplay of different interactions between atoms in disordered state. The values of A depend 
on the sign of pressure. 

PACS numbers: 45.05.+X, 61.43.-j, 62.50.-p 

Keywords: Disordered solids; Discrete systems; High-pressure effects 



I. INTRODUCTION 

As it was ascertained in last years a number of anoma- 
lous properties of glasses at low temperatures can be ex- 
plained supposing that there exist anomalous low-energy 
excitations: two-level (tunneling) systems (TLSs) p], [2[ 
and quasi-localized modes [3-5]. A characteristic prop- 
erty of these excitations is their remarkable sensitivity to 
different impacts including external pressure P. 

The pressure dependence of TLSs in glasses, on the 
basis of rather general arguments, was considered by 
Phillips 6] . He concluded that the anomalous (negative) 
value of the Griineisen constant for vitreous silica may 
be explained supposing that the TLSs, when pressure is 
applied, have a slight tendency to become less asymmet- 
ric. Such a tendency means that the energy difference 
between the two potential energy minima of a TLS A 
and the derivative of A with respect to P are to some 
extent correlated, so that the mean value (AdA/dP) 
differs from zero (and is negative). 

In 0], it was shown that the temperature cycling ef- 
fect on the width of spectral holes, burned in a dye-doped 
polymer glass at high pressure and low temperature can 
be understood if one assumes that under external pres- 
sure the higher minimum in the potential energy of TLSs 
gains more energy than the lower minimum. This means 
that in polymer glasses there also exists a correlation be- 
tween A and dA/dP; however, it is of the opposite (i.e. 
positive) sign, so that TLSs in these glassy materials be- 
come more asymmetric with increasing pressure. 

Although the existence of correlation between A and 
dA/dP in amorphous solids is expected, it is of interest 
to verify this property of disordered state by direct molec- 
ular dynamics simulation. For checking this correlation 
numerically, one must bear in mind that the concentra- 
tion of TLSs in amorphous state is usually rather low 
and therefore quite large clusters of atoms should be in- 
volved in computations, making them fairly intricate and 



time consuming. A relatively simple opportunity here is 
to compute one-dimensional (ID) disordered system con- 
sisting of a reasonably moderate number of atoms. 

In this study, we have carried out the computations 
of a disordered atomic chain. To produce a disordered 
state, special pair potentials having more than one en- 
ergy minimum were employed. In this model, already in 
a single-component case, one gets the states of the chain 
with interatomic bonds of different length. In our com- 
putations we employed the pair potentials of two types: 
the potential Zl in Ref. [8] with Friedel oscillations and 
the Schilling piecewise parabolic double- well potential [|| , 
which has a discontinuity of its first derivative at the bar- 
rier separating the minima. 



II. DISORDERED CHAIN WITH SMOOTH 
MULTI-WELL PAIR-POTENTIAL 

Let us first discuss a monatomic chain with the pair 
potential Zl [8j 



V(r) = a e ar cos (2Ax) /r 3 + b (a/r) n + V , 



(1) 
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where r is the distance between atoms in proper units. 
Values of the parameters a = 1.58, a = —0.22, b = 
4.2 ■ 10 s , a — 0.331, n — 18 are the same as those given 
in [8j for the potential Zl. The quantity Vq is a constant 
(we take Vo = 0). This potential has the sequence of 
minima and maxima describing the Friedel oscillations. 
This-type potentials are typical for metals where the os- 
cillations describe the effect of screening of the electric 
field by free electrons. The parameter k determining the 
period of oscillations is given by the Fermi wave vector. 
This parameter will be varied below to show that one 
can observe distinct kinds of response of the chain to the 
applied pressure. 

In our computations we examined the chain consisting 
of 400 atoms. We took into account the interactions be- 
tween 24 nearest atoms (inclusion of more interactions 
did not noticeably change the results). The disordered 
state was generated by randomly selecting 20 percent of 
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the nearest neighbor pairs of atoms and placing them at 
the distance, corresponding to the second minimum (sit- 
uated at r m 1.8) of the potential given by Eq. (1). All 
the other nearest neighbor pairs were initially placed at 
the distance (r ks 1.1), corresponding to the first (main) 
minimum. The initial configuration of atoms was allowed 
to relax to the local potential minimum of the chain. In 
order to find the relaxed configuration, the equations of 
motion of atoms were integrated using the forth-order 
Runge-Kutta algorithm. The motion towards the relaxed 
configuration was accomplished step by step, zeroing of 
the velocities of all atoms at every five to ten time steps. 
This procedure was repeated until only negligibly small 
changes (less than 10 -8 in our units) in the positions of 
all atoms were observed, indicating that the final config- 
uration of the chain was reached. 

TLSs in the present model are located at the bound- 
aries of "dense" and "rarefied" islands (groups of atoms 
with small (r k, 1.1) and large (r ss 1.8) distance, respec- 
tively). The tunneling transition in a TLS corresponds 
to the motion of the atom(s) located at the boundary 
of short and long bonds, resulting in the permutation of 
two bonds, i.e. in the displacement of a long bond to 
another position; see, e.g. the changes in the fragment 
of the atomic chain in the picture below - displacements 
of atoms (and bonds) cause the transitions between the 
configurations presented in the different lines of the pic- 
ture. 
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A. Free ends 

Our first computation was done for the chain with free 
ends. In the first run, the positions of atoms of the re- 
laxed initial configuration were determined and the po- 
tential energy of the chain and its length L$ were found. 
After the first run, one of the long bonds was displaced 
to one of the nearest positions (this displacement corre- 
sponds to the single-particle tunneling transition). Then 
the new coordinates of atoms in the relaxed configuration 
and the new length were computed and the new potential 
energy of the chain found. The difference of this energy 
value from the value corresponding to the initial relaxed 
configuration gave us the magnitude of A. The described 
runs were repeated with displacements of all other long 
bonds to the nearest positions, and the values of A of all 
other single-particle TLSs were found. 

Our computations revealed that the largest | A| belongs 
to the TLSs which correspond to the configurations pre- 
sented in the first and second lines of the picture. If we 
mark the short bond by 1 and the long bond by 2, then 
the back and forth tunneling transitions in these TLSs 



can be denoted as (1 1 2 2 1) o (1 2 1 2 1). The next 
largest |A| belongs to the TLSs which correspond to the 
configurations presented in the second and third lines of 
the picture (the tunneling transition in these TLSs can 
be denoted as (1 2 1 2 1) o (2 1 1 2 1)). The configura- 
tions with more than two adjacent long bonds, e.g., 
(1 2 2 2 1), (1 2 2 2 2 1), etc., lead also to the TLSs with 
a comparable |A|. The concentration of these structures 
in our case is the smaller the larger is the "rarefied" is- 
land. In the model under consideration, the structures 
with extra long bonds corresponding to the third, the 
fourth, etc. minimum of the potential given by Eq. (TTJ) 
can also exist. However, these structures were ignored 
relying on the physical argument of their weak stability. 

The energy difference A of a TLS depends also on 
the degree of disorder. This dependence is caused by 
the interaction of the described above central atoms of 
the TLS with the surrounding atoms. This conclusion 
is supported by consideration of the chains with lower 
and higher concentration of long bonds. Namely, it was 
found that in the first case the dispersion of the values of 
A of the TLSs with the same configuration of the central 
atoms is diminished and in the second case enlarged. 

We have also studied the cooperative TLSs that corre- 
spond to a simultaneous displacement of pairs and larger 
groups of atoms (e.g. tunneling transitions (112 2 11) 
<-> (2 1 1 2 1 1), (1 2 1 2 1 2) o (2 1 2 1 2 1), etc). It was 
found that the main conclusions which were drawn above 
for single-particle TLSs, hold also for cooperative TLSs. 
However, as more particles have to move, the effective 
mass of the tunneling group of atoms is increased, which 
implies rather long relaxation times. Therefore, the con- 
tribution of such TLSs to low-temperature characteristics 
can be neglected on a time scale much shorter than the 
corresponding relaxation time f]~oj ] . 



B. Isobaric and isochoric tunneling transitions 

The difference of the lengths of different configurations 
of the chain with free ends is very small as compared to 
Lq. The larger the Lq is the smaller is the pressure and 
the work required for bringing the lengths of these con- 
figurations to the same length Lq. Thus, for sufficiently 
large Lq, the energies of the configurations of the chain 
with free ends are practically equal to the energies of 
the corresponding configurations of the chain with fixed 
length Lq. The same is true for the corresponding values 
of A. 

An analogous situation exists for isobaric and isochoric 
tunneling transitions for nonzero pressure. An external 
pressure applied to the ends of the initial configuration 
of the chain causes the change of the length of the chain 
from Lq to L. In the isobaric case (the pressure P is 
fixed) after a tunneling transition the latter length in its 
turn changes from L to L + dL. The difference A of the 
energy of corresponding configurations includes the term 
PdL. The work required to bring these configurations to 
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the same length L, i.e. to the isochoric case, is also PdL 
if we assume the change of P required for bringing the 
lengths of these configurations to the same length L to be 
negligibly small. The larger the chain is the better this 
assumption is fulfilled, and the smaller is the difference 
between the isobaric and isochoric A values. In our case 
of the chain of 400 atoms the difference is less than 0.1 
%. 

Here, following 6] we consider the isochoric tunneling 
transitions. We take pressure P in reduced dimensionless 
units according to the relation P = (L — L)/L . The 
case P — corresponds to the free ends of the initial 
(i.e. prior to the tunneling transition) configuration of 
the chain. 



C. Effect of compression 

In the next step the effect of positive pressure (in our 
case a uniaxial compression of the chain) was examined. 
Here the positions of the end atoms were fixed so that the 
length of the chain L would be the same (in the P = 
case) or shorter (in P > case) than the length Lq of the 
initial relaxed configuration of the chain with free ends. 
The initial positions of all other atoms were chosen so 
that the distances between atoms would be uniformly 
reduced. By using this initial condition, the relaxed po- 
sitions of all by end atoms were computed and the poten- 
tial energy of the chain found. Then in a similar manner 
as in the runs without pressure, one of the long bonds 
was displaced to another (nearest) position, the new re- 
laxed configuration was computed and the new potential 
energy found. The difference between this energy value 
and the value corresponding to the initial relaxed config- 
uration gave us the energy difference of the two minima 
of a TLS at a pressure P, denoted by A(P). The de- 
scribed runs were also repeated with displacements of all 
other long bonds, and the values of A(P) of all other 
TLSs were found. 

Figures 1 and 2 show the obtained dependence of 
A(P) - A(0) on A(0) for different TLSs at a small pres- 
sure P = 0.001 in the case of the potential Zl given by 
Eq. (JIJ with the parameters k = 4.12 and k = 4.5, re- 
spectively. From these figures it is clearly seen that A(0) 
and A(P) — A(0) are correlated. The sign of the corre- 
lation may be different: in the case k = 4.12 (Fig. 1) 
the correlation is positive (TLSs become more asymmet- 
ric under pressure), in the case k = 4.5 (Fig. 2) it is 
negative. The correlation is almost complete. However, 
there exist intermediate cases (with overall weak depen- 
dence of A on P) where the correlation is insignificant 
(see, e.g., Fig. 3). 

Note further that in the case k = 4.12 the minimum 
of a TLS with a higher energy corresponds to a larger 
length of the chain with free ends. However in the case 
k = 4.5 the relation between energy and length is just the 
opposite: the minimum with larger energy corresponds 
to a smaller length. In both cases the larger is the length, 
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FIG. 1: Changes in A of TLSs under dimensionless pressure 
P — 0.001 in the disordered atomic chain with the pair po- 
tential Zl. The presented case corresponds to k = 4.12 in Eq. 
(fTJ); the other parameters are given in the text. 
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FIG. 2: The same as in Fig. 1, but for k = 4.5. 



the bigger is the change (increasing) of the energy with 
pressure. This can be understood as follows - to compress 
the chain extra energy is needed, which is the larger the 
larger is compression. 

Unlike the case of low P values, in the case of high P 
values the correlation between A and OA/dP was found 
to be always positive. In the region of high enough pres- 
sures the TLSs are more asymmetric than at zero pres- 
sure. It is clearly seen in Fig. 4, where the dependence 
of the energy differences A of two TLSs of the types 
(1 1 2 2 1) (1 2 1 2 1) and (1 2 1 2 1) (2 1 1 2 1) on 
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FIG. 3: The same as in Fig. 1, but for k = 4.3. 



dimensionless pressure P is given for the chain with the 
same interaction parameters as in Figs. 1 and 2. The 
upper values of P in Fig. 4 correspond to the critical 
pressure at which the abrupt changes in a chain begin to 
occur. Such changes take place when some long bonds 
transform into the short bonds, which results in disap- 
pearance of respective TLSs. Obviously, this happens in 
such a manner that TLSs gradually increase their asym- 
metry. With further increase of pressure more and more 
TLSs disappear. At sufficiently high pressure P ~ 0.1 
there remain no TLSs at all and the chain becomes or- 
dered. 

The presented results are in agreement with the ex- 
periments 

GH3, 

where a reduction in the number of 
TLSs and soft localized modes with pressure in glasses 
was observed. The reduction in the number of states of 
a disordered structure with the growth of its density was 
earlier found theoretically in [14j by an analytical study 
of a chain of particles with Schilling piecewise parabolic 
double-well potential An analogous result was ob- 
tained in [15] by the molecular dynamics simulation of 
32 and 108 particles with the finite-range Lennard- Jones 
interaction between them. 



D. Effect of stretching 

It is reasonable to expect that the response of a glass to 
the volume contraction and dilatation should be different 
[l6l Il7j |. If so then the value of A(P) should depend on 
the sign of P. In order to control this statement we 
performed computations of TLSs also for negative P (in 
our case for uniaxial stretching of the chain). We have 
found that, indeed, the value of the energy A of a TLS 
depends on the sign of pressure, which can clearly be seen 
in Fig. 4. Only in the case of a large \P\ the change 
of A with \P\ for positive and negative P is similar: 
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FIG. 4: Pressure dependence of |A| in the disordered atomic 
chain for the two TLSs with the largest |A| values. Pair po- 
tentials are the same as in Figs. 1 and 2; solid lines correspond 
to k — 4.12 while dashed lines to k = 4.5. 



in both cases A increases with \P\. It was also found that 
at negative pressure P w —0.02 the chain breaks down. 



III. DISORDERED CHAIN WITH PIECEWISE 
PARABOLIC DOUBLE- WELL POTENTIAL 

Lastly we have performed computations of the pressure 
dependence of ID disordered chain with the Schilling- 
type piecewiseparabolic double- well potential for nearest 
neighbors d, Hj] 

Vi(r)= K,^i(r-r 1 ) 2 e(r -r) + [ a2 (r-r 2 ) 2 +6]O(r-r ), 

where 0(x) is the Heaviside step function. In [t| [l^ 
only the case a\ — a 2 was considered. Here the values 
of parameters Vq = —0.743, a\ — 40, a 2 = 5.4, b = 0.6, 
?"o = 1-357, t*i = 1.13, r 2 = 1.88 were chosen such that 
the minima of the potential given by Eq. ([2]) would coin- 
cide with the first two minima of the used above potential 
Zl. In the model [t| [lij], besides the interaction between 
nearest neighbors, also the interaction with next-nearest 
neighbors is taken into account. This interaction is rep- 
resented by a parabolic potential 



V 2 (r) = c(r-r 3 ) 2 



(3) 



We took c = 0.05 or c~ —0.05 and r 3 = 4.0. The results 
of computations of TLSs for small dimensionless pressure 
P = 0.001 are shown in Fig. 5. It is clearly seen that 
the quantities A(0) and A(P) — A(0) are correlated, and 
the correlation is complete and negative. The variance of 
the TLSs parameters is very small - much smaller than 
in the case of the potential Zl. This is a consequence of 
a strong localization of the Schilling interactions. 
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FIG. 5: Changes in A under dimensionless pressure P — 0.001 
in a disordered chain with the Schilling-type pair potential. 
The presented cases correspond to c = 0.05 and c = —0.05 in 
Eq.©; the other parameters are given in the text. 



We have also performed computations for other values 
of the parameters of the potentials Vi(r) and V-^r). It 
appeared that the sign of the difference a2—a\ determines 



the sign of the correlation between A and dA/dP: ifai > 
ci2 then the sign of the correlation is negative, thereby if 
a\ < a2 then the sign is positive. In the case a\ = ai the 
values of A do not depend on P. In the large P limit all 
the TLSs disappear. 



IV. CONCLUSION 

Molecular dynamics simulations of ID disordered state 
unambiguously testify that for two-level systems there 
exists a correlation between the magnitudes of A and 
their alterations under external pressure. The correla- 
tion may be either positive or negative sign. This sig- 
nifies that the asymmetry of two-level systems may in 
average increase or decrease with pressure, depending on 
the interplay of interactions between particles in the dis- 
ordered state. The response of TLSs to external pressure 
depends on the sign of pressure. 
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